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Abstract 



At zero temperature magnetic phases of the quantum spin-1/2 Heisenberg antiferromagnet on a 
simple cubic lattice with competing first and second neighbor exchanges (Ji and J2) is investigated 
using the non- linear spin wave theory. We find existence of two phases: a two sublattice Neel 
phase for small J2 (AF), and a collinear antiferromagnetic phase at large J2 (CAF). We obtain the 
sublattice magnetizations and ground state energies for the two phases and find that there exists 
a first order phase transition from the AF-phase to the CAF-phase at the critical transition point, 
Pc = 0.28. Our results for the value of pc are in excellent agreement with results from Monte-Carlo 
simulations and variational spin wave theory. We also show that the quartic 1/5" corrections due 
spin-wave interactions enhance the sublattice magnetization in both the phases which causes the 
intermediate paramagnetic phase predicted from linear spin wave theory to disappear. 




PACS numbers: 75.10.Jm, 75.40.Mg, 75.50.Ee, 73.43.Nq 



I. INTRODUCTION 



Frustrated quantum Heisenberg magnets with competing nearest neighbor (NN) and 
next-nearest-neighbor (NNN) antiferromagnet (AF) exchange interactions, Ji and J2 re- 
spectively, have been under intense investigation both theoretically and experimentally in 
condensed matter physics for more than a decade.- At low temperatures these systems ex- 
hibit new types of magnetic order and novel quantum phases.-i^>^ A well-known example is 
the quantum spin-1/2 antiferromagnetic Ji — J2 model on a square lattice, which has been 
studied extensively by various analytical and numerical methods 1^^^'^'^^^^^^ For this 
two-dimensional square lattice system with J2 = the ground state is antiferromagnetically 
ordered at zero temperature. Addition of next nearest neighbor interactions induces a strong 
frustration and break the antiferromagnetic (AF) order. The competition between the NN 
and NNN interactions for the square lattice is characterized by the frustration parameter 
p = J2/ Ji- It has been found that a disordered quantum spin liquid phase exists between 
Pic ~ 0.38 and p2c ~ 0.60. For p < pic the square lattice is AF-ordered whereas for p > p2c 
a coUinear phase emerges. In the collinear state the NN spins have a parallel orientation in 
the vertical direction and antiparallel orientation in the horizontal direction or vice versa. 
The nature of phase transition from AF-ordered state to disordered state at pic is of second 
order and from the disordered state to the collinear state at p2c is of first order. 

The properties of quantum magnets depend strongly on the lattice dimensionality since 
the tendency to order is more pronounced in three dimensional (3D) systems than in the 
lower dimensional systems. Furthermore, in 3D the available phase space is more and we 
expect quantum fiuctuations to play a lesser role as compared to ID and 2D. In ID and 
2D the available phase space is limited and quantum fiuctuations play a dominant role in 
determining the quantum critical points. Despite this fact a magnetically disordered phase 
has been observed in frustrated 3D systems such as the Heisenberg AF on the pyrochlore 
lattice^^ or on the stacked kagome lattice^»^>^»^. Studies on the Heisenberg AF on the 
pyrochlore lattice (a geometrically frustrated system) have revealed the existence of a spin 
liquid state.~ On the other hand, for the 3D J1-J2 model on the body-centered cubic 
(BCC) lattice there are no signs of an intermediate quantum paramagnetic phase at zero 
temperature.— For the BCC lattice competing interactions and not lattice geometry 
generates the frustration. This comparison illustrates how the magnetic phase diagram may 
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dramatically change based on whether the frustration is generated by competing interactions 
or by geometry. 

Most of the efforts on quantum 3D magnets have primarily focused on geometrically 
frustrated lattices.-"^ There exists some computational2^>^'2^>2ii2^i2^ZLiS^;2Ei^;^i2^!2^ and very 
few analytical studies^^i^SiMiS^ of the magnetic phase diagrams and magnetic order of spin- 
1/2 Heisenberg AF on 3D lattices where on the study of magnetic phase diagrams and 
magnetic order of spin-1/2 Heisenberg AF on 3D lattices where competing interactions induce 
frustrationM-^^-^ '^^'^^'^^^'^'^' '^^''^^'^^'^^'^--^^-^ 

Very few analytical and numerical results exist for the the frustrated Ji — J2 isotropic 
Heisenberg model on a simple cubic lattice.— >^i^>^ This model has been studied previously 
using Monte Carlo simulations^, variational spin wave theory^^, and modified spin wave 
theorj*i^. In a recent work the critical properties of the 3D anisotropic quantum spin-1/2 
model on a simple cubic (SC) lattice has been investigated within the framework of the 
differential operator technique and by using an effective field theory in a two-spin cluster.— 
The study revealed that at zero temperature there is a AF-lamellar (first order) phase 
transition. The motivation for the present work is to investigate the zero temperature 
phases of this model in the framework of non-linear spin wave theory (NLSWT) and to 
obtain the critical transition points of this model. Also we will compare our results from 
NLSWT with the prediction from the linear spin wave theory (LSWT). 

The paper is organized as follows. In Section[niwe set-up the Hamiltonian for the spin-1/2 
Heisenberg AF on the SC lattice. The classical ground state configurations of the model and 
the different phases are then discussed. In Section IIIII we map the spin Hamiltonian to the 
Hamiltonian of interacting bosons and develop the NLSWT sublattice magnetization and 
energy expressions. The sublattice magnetizations and the ground state energies for the two 
phases are numerically calculated and the results are plotted and discussed in Section HVl 
Finally we summarize our findings in Section |Vl 
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II. CLASSICAL GROUND STATE CONFIGURATIONS 



The Hamiltonian for a spin-1/2 Heisenberg antiferromagnet with first and second neighbor 
interactions on a simple cubic (SC) lattice is 

h = Ij^Y.s,-s, + ^j,Y,s,-s,. (1) 

(«J> [ij] 

Ji is the NN and J2 is the frustrating NNN (which are along the face diagonals of the cube) 
exchange constants. Both couplings are considered antiferromagnetic, i.e. Ji, J2 > 0. For 
the SC lattice the number of nearest and next-nearest neighbors are Zi = 6 and Z2 = 12. 

The limit of infinite spin, S —>■ 00, corresponds to the classical Heisenberg model. We 
assume that classically the spin configurations of the system are described by Si = 
where u is a vector expressed in terms of an arbitrary orthonormal basis and q defines the 
relative orientation of the spins on the lattice. The classical ground state energy of the 
lattice in terms of the frustration parameter, p, is given by 

^k/iV^i = ^^'^i[7ik + P72k], (2) 



with the structure factors 



7ik = 3 
1 

72k 



cos(/cx) + cos{ky) + cos{kz) , (3) 
cos{kx) cos{ky) + cos(A;j^) cos(A;2) + cos(A;2) cos{kx) , (4) 
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where we define the parameter of frustration as p = Z2J2/ ziJi-— The wave- vectors along the 
X, y, and z directions are denoted by k^, ky, and kz. The number of lattice sites are given 
by and we have set the lattice spacing a = 1. 

At zero temperature the classical ground state (GS) for the SC lattice can be characterized 
by the values of p. p = corresponds to the unfrustrated case (only AF interactions between 
NN) . For p < 1/2 or J2/ Ji < 1/4, there is a single minimum in energy Eq / NJi = —33"^ ( 1 —p) 
for the wave- vector (±7r, ivr, ±7r). They correspond to the classical two sublattice Neel state 
(AF phase) where spins in A and B-sublattices point in opposite directions [Fig. l(a)| . 

For p > 1/2 apart from the global rotation the classical ground state has an infinite degen- 
eracy - the frustration is uniformly distributed on all the spins, causing a non-coUinear GS 
with very large degeneracy. In general, the GS for p > 1/2 can be decomposed into two NNN 
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tetrahedra. The spin configurations in each of these two NNN tetrahedra can be character- 
ized by two angles 9 and (f). This results in a four sublattice A, B, C, D [Fig. l(b)| antiferro- 
magnetic structure. Out of these infinite possibihties there are three collinear configurations 
(one hne up, one hne down). The wave- vectors corresponding to these colhnear states are 
(7r,7r,0), (7r,0,7r), (0,7r,7r). The classical GS energy for these states is -Eq/A^Ji = —5*^(1 
Thermal or quantum fluctuations lift these degeneracies and select specific discrete states 
and it has been conjectured that thermal or quantum disorder favors collinear states {order 
by disorder) .^^'"^"^ The four fold rotational symmetry of the lattice is spontaneously broken 
in this state. By employing a spin wave theory based on the general four sublattice mean 
field ground state it has been shown that the quantum fluctuations stabilize a collinear spin 
ordering.-^ Quantum Monte Carlo simulations on the frustrated SC lattice for p > 1/2 also 
confirm this conjecture.— In the present article, for p > 1/2, we consider the system to be 
in one of these three collinear configurations (collinear antiferromagnet or CAF). 

p = 1/2 corresponds to the case where both Ji and J2 compete - causing frustration in 
the system. This critical value Pciass = 0.5 is the classical phase transition point where a 
phase transition from AF to CAF phase occurs. In this work we will investigate the role 
of quantum fluctuations in the two different phases (AF and CAF) of the model and how 
these fluctuations shift the critical transition point. 



III. SELF-CONSISTENT NON-LINEAR SPIN WAVE THEORY 

The Hamiltonian in Eq. [1] can be mapped into an equivalent Hamiltonian of interacting 
bosons by transforming the spin operators to bosonic operators a, and b, b^ using the well- 
known Holstein-Primakoff transformations. For the AF-phase ( J2 < Ji) the operators a, 
and b, W are for the A and B sublattices. On the other hand for the CAF phase ( J2 > Ji) 
we have used the operators a, and 6, 6^ for the up and down spin configurations. 

- V^(l-^)a., S-^,^V2S4[l-''^), 

S j^j^ — — S dj^di^ 

S^^ ^ V2Sb](l - ^), S,^ ^ V2S(1 - 

SI,, = -S + b]b„ (5) 
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In these transformations we only kept terms up to the order oi 1/S. Next using the Fourier 
transforms 

k k 

the real space Hamiltonian is transformed to the k-space Hamiltonian. In the following two 
sections we study the cases J2 < Ji and J2 > Ji separately. 



1. J2 < Ji- A F phase 



In this phase the classical ground state is the two-sublattice Neel state [Fig. l(a)| . For 



the NN interaction, spins in A sublattice interacts with spins in B sublattice and vice versa. 
On the other hand the NNN exchange J2 connects spins on the same sublattice, A with A 
and B with B. Substituting equations ([5]) into ([T]), the k-space Hamiltonian takes the form: 

if = i^W + + ^(4) (6) 

The classical ground state energy H^^^ and the quadratic terms H^"^^ are 

= -InJ,S'z,{1~p), (7) 
= J^Sz^ J2 [4'k (4«k + blb^) + 5jLVjc&-k + ^-M] ■ (8) 



with the coefficients A^^^ and B^y^ defined as 

4k = l-p(l-72k), (9) 
< = 7ik. (10) 

The quartic terms in the Hamiltonian H^'^^ involve interactions between A — B (for NN 
terms) and A — A, B — B (for NNN terms) sublattices. The Hamiltonian for these interaction 
are stated in Appendix \^ Eq. IA1[ These terms are evaluated by applying the Hartree-Fock 
decoupling process. The contributions of the decoupled quartic terms to the harmonic 
Hamiltonian in Eq. [8] are to redefine the values of A^^^ and B^ which are now 

^« = (i-!fi±^)-p|i-„,|(i-!fi^), (11) 



0.(1-^), (12, 



5k ...^ s 



^k 



(1) _ .lfAW\^ fnW^^ 



A\^> - Bl^> . (13) 



The coefficients ui,vi and wi are in Appendix |X1 They are evaluated self-consistently from 
equations (HH), (HI, flM|) - fE4l) . 



The quartic corrections to the ground state energy is calculated from the four-boson 
averages. In the leading order they are decoupled into the bilinear combinations (equations 
(IA2p - (lA4p ) using Wick's theorem. The corresponding four boson terms are, 



(alaib^jbj) = ul + v'f, {aJMhjhj) = 2uiVi 



{ajaiaibj) = 1u\V\, {alaiOjaj) = uf + w'f, 
{aia}ja}-aj) = 2uiWi, (ajaiaia^j) = 2uiWi. 

This yields the ground state energy correction from the quartic terms: 

5^(4) = -InJiZi [{ui + - p{ui - w,f] . 

Adding all the corrections together the ground state energy takes the form 



E/NJi 



-\z,S{S+l){l-p) + \z,S 



1 



N 



E 



(1) 



+ [("1 + '"^^^^ - Ml - Vi) - p{ui - Wi){l - Ml + Wi)] 



and the average sublattice magnetization (So) is given by 



{So) — s 



-{- 

2S\n 



A 



(1) 



\^ Ilk 



k ^k 



(1) 



(14) 



(15) 



(16) 



(17) 



Using equations (fTTl) - (IT^ . we numerically evaluate E/NJi and {So}- 



2. J2 > Ji .• CAF phase 

The classical ground state for J2 > Ji is considered to be in one of the three coUinear 
states [Fig. |l(b)| . For NN and NNN exchanges there are A - B, A - C , A - D, B - C , 



B — D, and C — D interactions between the four sublattices [See Fig. l(b)| . Considering all 



their contributions together up to the quadratic terms the harmonic Hamiltonian takes the 
same form as before with 

= -\nJ,S^z^{1+p), (18) 
^ok = ^(1 + cos/cz) + ip(l + cos/ci^cos/cj^), (19) 
-^ok — -;:{coskx + cosky){l + pcoskz). (20) 



The quartic terms in the Hamiltonian for this case are shown in Appendix [Bl These terms 
are decoupled and evaluated in the same way as before. The renormalized values of the 
coefficients A^"* and B^^ are 

+ p{(u — wi){l — cos cos ky) — 2(ll + V2)}], (21) 
= ^ok ~ ^{("^ + ^^1) + ^^2) cosA;2}(cosA;^ + cosfcy), (22) 

The coefficients u, vi,V2,wi,Wz are in Appendix [Bl As before these coefficients are calculated 
self-consistently from equations fl^ - (|25]) and (]B2p - (]B8p . The quartic correction to the 
ground state energy (following the same Hartree-Fock decoupling process as done in the 
AF-case) is 

5^(4) = -InJiZi \2(u + viY - (m - + 2p(u + V2Y - piu - W2Y] . (24) 

D 

Combining all these corrections, the ground state energy takes the following form: 
E/NJ, = -y,S{S+l){l+p)+^-z,s[^J2^^^^ 

k 

- ^zi [(u - w,) - 2(u + vi) + p(u - wi) - 2p(u + V2)] 

- ^z, [2{u + vif -(u- + 2p(u + V2f - p(u - W2f] . (25) 

The sublattice magnetization and the ground state energy are then obtained numerically 
using equations f|T71) and fl25|l . 



IV. RESULTS 

In Fig. [2] we show the self-consistent values of the different parameters Ui,Vi,Wi (AF 
phase) and Ui,Vi,V2,Wi,Wz (CAF phase) of our model. These parameters which provide the 
quartic corrections to our model do not appear in the LSWT calculations for the sublattice 
magnetization, {Sa) and the ground state energy, E. We see from Fig. [2] that most of these 
coefficients vary significantly with p especially as p approaches 0.5 from both ends. This 
demonstrates that non-linear corrections due to the spin-wave interactions play a significant 
role in determining the different phases of our model. 
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Figure [3] shows the result for the average sublattice magnetization for the SC lattice for 
both AF and CAF phase without (dashed line) and with (solid line) quartic corrections. In 
the AF ordered phase or the two sublattice Neel phase where A and B sublattice spins point 
in the opposite directions, sublattice magnetization decreases monotonically with increase 
in p until p ~ 0.49. This gradual decrease in (So) with increase in p is expected as increasing 
strength of NNN interaction J2 disorders the antiferromagnetic spin alignments. With only 
quadratic terms in the Hamiltonian (linear spin wave theory) we find that (Sa) approaches 
zero as p Pd where Pd ~ 0.48 indicating a order-disorder phase transition to a disordered 
paramagnetic (PM) state at this point. In the CAF phase with lines of spins up and down, 
LSWT calculations show that that {Sa) decreases as p approaches 0.5 from above and at 
P = Pc2 = 0.50 there is an another phase transition from the CAF state to the disordered 
PM state. This is similar to the two dimensional AF-square lattice with Heisenberg spins 
where we have a line of quantum critical points between 0.38 < p < 0.60. However, self- 
consistent calculations with quartic 1/S corrections drastically alter the zero temperature 
phase diagram. We find that in the AF-phase with increase in p the system aligns the 
spins antiferromagnetically along the horizontal and vertical directions - thus decreasing the 
sublattice magnetization from ^ 0.42 for p = to ^ 0.30 for p = 0.49. In the CAF phase 
(So) steadily decreases from ^ 0.41 for p = 1 to ~ 0.27 for p = 0.52. There is no existence 
of any disordered state as predicted by the linear spin- wave theory (quadratic corrections). 
The disordered PM region disappears completely and we only obtain two phases: AF and 
CAF. This is one of our main findings in the present work. This significant change due to 
the quartic corrections is due to the enhancement of order by quantum fluctuations. 

At p = (no frustration) there is no quartic corrections to {Sa)- This can be observed 
from equations ffTTl) - f|T3l) as the correction factor (1 — {ui + vi)/S) cancels out in equation 
[T71 Our non-linear spin wave theory calculations become unstable close to the classical 
transition point Pciass = 0.5 since the coefficient A^^ becomes equal to S^- 

We have also applied the NLSWT technique to compute the quartic corrections in the 
spin-1/2 Heisenberg AF on a body-centered lattice.— The LSWT calculation for the BCC 
lattice does not predict any intermediate disordered state and the quartic corrections play a 
role in stabilizing the sublattice magnetization (see Fig. 2 of Ref. 1221). However, the effect of 
quartic corrections is more pronounced in the SC lattice where the intermediate disordered 
phase disappears. 
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In Fig. m we show the ground state energy per site, E/NJi, for the AF and the CAF 
phases with (sohd hne) and without (dashed hne) quartic corrections as a function of the 
frustration parameter p. Classically Pciass = 0.5 or J2/J1 = 0.25 is the critical point where a 
phase transition from the AF phase to one of the three CAF phases occur. With increase in 
frustration (as p approaches 0.5 from both sides) we expect the GS energy to increase as the 
system goes from an energetically favored ordered state to a more disordered state. However, 
linear spin-wave theory calculation fails to capture this. Especially when p is close to 0.5 we 



find a slight downward turn in energy. This has been reported in Ref. |38|. On the other hand. 



NLSWT calculation correctly produces the expected energy increase. At p = the calculated 
energy with the quartic corrections is slightly lower than the energy obtained without the 
quartic corrections. This small decrease from the LSWT calculation is due to the ground 
state energy correction, which is negative (as seen in Eq. [15] - these terms originate from 
the self-energy Hartree diagrams). As our spin- wave theory calculation becomes unstable 
in the regime 0.49 < p < 0.52 we used a spline fit for the AF-phase energy data points and 
then extrapolated the line so that it intersects the CAF-phase energy line. The extrapolated 
curve is shown by dotted lines (color online) in the figure. After extrapolation, we find that 
the two energies meet at pc ~ 0.56 or J2/J1 ~ 0.28. The symmetries of the two phases are 
different: SO(3)/SO(2) for the AF phase and Z3xSO(3)/SO(2) for the CAF phase.^ Due to 
the different symmetries of the two phases the transition is of first order. This is confirmed 
by the kink in the energy at pc ~ 0.28. Our obtained value from our self-consistent NLSWT 
calculations for the quantum critical point is J2/ Ji ~ 0.28. This is another major finding of 
our work. 



Using the variational spin- wave theory the authors in Ref. |38| obtained an upper bound of 
0.27 for the ground state energy. However, their variational calculation slightly overestimates 
the value of GS energy. This is quite noticeable at p = 0, where their energy value is higher 
than the LSWT prediction. We have found the p = energy value to be slightly less than 
the LSWT prediction. This is due to the quartic corrections explained earlier. The authors 
in Ref. Q used a modified spin-wave theory based on Dyson-Maleev representation of the 
spins to study this model. They numerically obtained the value of critical transition point 
to be Pc = 0.30. The other known existing numerical work is by Diep et. al.— By extensive 
standard and histogram Monte-Carlo simulations, they obtained the transition point to be 
around 0.26. Our result is in excellent agreement with the results obtained from Monte- 
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Carlo and variational spin-wave theory calculations (our pc differs by less than 3.5% from 
the variational spin- wave theory prediction). 

V. CONCLUSION 

In this work we have investigated the zero temperature phases of a spin-1/2 Heisenberg 
frustrated AF on a SC lattice by considering the quartic 1/ S corrections due to the spin- 
wave interactions. We have compared our results obtained from NLSWT calculations with 
the predictions from LSWT. It is known that LSWT predicts the existence of three phases: 
a two sublattice Neel phase for smaller values of the NNN exchange J2, an intermediate 
paramagnetic phase, and a coUinear phase for larger values of J2. At zero temperature there 
are two quantum phase transitions - one from the AF-statc to the disordered paramagnetic 
state and the other from the disordered state to one of the three collinear states. Both these 
transitions occur at the quantum critical point Pc ~ 0.5 or J2/J1 ~ 0.25. We have found 
that quartic corrections significantly alter this phase diagram as intermediate paramagnetic 
phase disappears. We find the existence of two phases at zero temperature: a two sublattice 
AF phase for small J2 and a collinear phase phase for large J2. With the inclusion of 
quartic interactions the intermediate paramagnetic phase completely disappears. Due to 
the different symmetries of the two phases (AF and CAF) the transition between the two 
phases is of first order and we find that the critical point for transition to be 0.28. Our 
obtained result is in excellent agreement with existing numerical results from Monte Carlo 
simulations. 
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APPENDIX A: QUARTIC TERMS FOR THE AF PHASE 

The quartic terms for the AF phase from the NN interactions involve A — B interactions 
and for the NNN interactions involve A — A and B — B type interactions. Considering all 
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these interactions the quartic Hamiltonian takes the form: 



H^'^^ = - Jl X ["ajajfejfej + - (aibjbjbj + a\aiaihj + h.c 



{ij) 



+ -J2 X] |^a]^aja^aj ^Ojajajaj + alaittiOj + h.c^ + a ^ b . (Al) 

[ij] 

In the harmonic approximation the following Hartree-Fock averages are non-zero for the 
SC-lattice Heisenberg antiferromagnet: 



where ui. 



(1) 



1 r 2 



u 



1 = {4a.) = m)=-^[^j:^-^ 



A 



(1) 



k ^k 



vi = {aibj) = {alb]) 



ir2_^7ik5^ 



k ^k 



1 r 2 



2 LA^ 



E 



72k^: 



(1) 



k ^k 



(1) 



A 



(1) 



5, 



(1) 



(A2) 
(A3) 
(A4) 



APPENDIX B: QUARTIC TERMS FOR THE CAF PHASE 

The quartic terms for the collinear phase from the NN interactions involve interactions 
between the sublattices B — C and A — D along the x axis, A — C and B — D along the y axis, 
and A — B and C — D along the z axis. For the NNN interactions spin-spin interactions are 
between sublattices A — B and C — D (in the x?/-plane), A — C and B — D (in the a:2;-plane). 



and A — D and B — C (in the yz-^Aane) as shown in Fig. 1(b) Adding all the contributions 
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together yield 

{ijy} 



(ijyz) 



J, 



(ijxz) 



(Bl) 



Above jx,jy,jz are NN lattice sites along x,y,z axes and jxy,jyz,jxz connects one lattice 
site with a NNN corner lattice sites on the xy, yz, xz planes. The different coefficients that 
originate from Hartree-Fock decoupling process are 



(2) 



k 



(2) 



- 1 



Vlx = 



1 r 2 cos /cx-Bk 



(2) 



k '^k 



(2) 



1 r 2 ■^-^ cos kyB^ 



(2) 



k ^k 



(2) 



V2yz = (ai&j,.) = 



1 r 2 cos ky cos /c^-B,^ 



1 r 2 cos kx cos kzB^ 



E 

k 

(2) 

^(5) J' 



(2)-, 



(2) 



(2) 



(2) 



1 r 2 ^ cos kzAy^ 



.(2) 

k "^k 



1 r 2 cos kx cos kyA^ 



(2) 



a;. 



(2) 



(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
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where u;^^'' = J (^k^) - (^k ^) • % symmetry Vi^ = Viy = Vi and V2yz = V2xz = V2. 
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FIG. 1: AF and collinear antiferromagnetic (CAF) ordered phases of the SC lattice. In the AF 
phase all ^-sublattice spins point in the direction of an arbitrary unit vector while B-sublattice 
spins point in the opposite direction. For the CAF phase the spin configurations (lines of spins 
up and down) of the four sublattices A,B,C, and D are shown in the Figure. There are two other 
equivalent configurations with lines along the two other directions of the cubic lattice. 
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FIG. 2: Self-consistent results for the different parameters, ui,vi,wi (for the AF-phasc) and 
ui,vi,V2,wi,Wz (for the CAF phase) are plotted with the frustration parameter p = Z2J2/Z1J1 
(for the SC lattice zi = 6 and Z2 = 12). These coefficients vary significantly with p, which shows 
that the quartic interaction terms play a significant role in determining the different phases of our 
model. 
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FIG. 3: Average sublattice magnetization, (Sa), is plotted with the frustration parameter p for AF 
and one of the three CAF phases with (sohd hnes) and without (dashed hues) quartic corrections. 
At zero temperature without the quartic 1/S' corrections (Unear spin- wave theory) (Sa) — >■ at 
Pel ~ 0.48 indicating a phase transition from the AF-ordered state to the disordered paramagnetic 
state. At Pc2 = 0.50 there is a second phase transition from the colhnear state to the disordered 
state for T = 0. Non-hnear spin wave theory provides significant corrections to this phase diagram. 
With the quartic 1/S corrections the disordered PM region disappears completely and we only 
obtain two phases: AF and CAF. There is no existence of any disordered state as predicted by the 
linear spin- wave theory (quadratic corrections). For both the phases the quartic corrections to the 
Hamiltonian enhance the magnetic order. 
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FIG. 4: Zero temperature ground state energy per site, E/NJi, is plotted versus p without (dashed 
hnes) and with (sohd hues) quartic corrections for both AF (p < 0.5) and CAF (p > 0.5) ordered 
phases. Spin wave theory becomes unstable close to the classical transition point (pciass ~ 0.5) 
between the two phases. After extrapolation (shown by the dotted line), we find that the two 
energies meet at the quantum transition point, pc ~ 0.56 or J2/J1 ~ 0.28. This kink in the energy 
indicates a first-order quantum phase transition from the AF to CAF phase. Compared to the 
results without quartic corrections (long dashed lines) we find that the quartic corrections provide 
significant corrections to the ground state energy especially near the AF-CAF phase transition 
point. 
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